Anomalous dynamics of confined water at low hydration. 
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The mobility of water molecules confined in a silica pore is studied by computer simulation in 
the low hydration regime, where most of the molecules reside close to the hydrophilic substrate. A 
layer analysis of the single particle dynamics of these molecules shows an anomalous diffusion with 
a sublinear behaviour at long time. This behaviour is strictly connected to the long time decay of 
the residence time distribution analogously to water at contact with proteins. 



I. 



INTRODUCTION 



Water plays a major role in many different biological, 
chemical and physical phenomenai^. In most of these 
cases a large fraction of water is at contact with different 
substrates and its motion is restricted in small spaces. 

It is expected that both the geometrical confinement 
and the interaction with the substrate perturb the struc- 
tural and dynamical properties of water and some gen- 
eral trends have been found in experiments and computer 
simulation^. 

In many phenomena, like those connected with biolog- 
ical matter, it is of great relevance the behaviour of the 
shells of water close to substrates. In particular the slow 
dynamics of water close to the surfaces of proteins might 
play a fundamental role in the protein functionality and 
evidences have been found of a glasslike behaviour of wa- 
ter at contact with plastocyanin investigated in a wide 
temperature range by computer simulation^. 

Water confined in Vycor glass is a good prototype sys- 
tem for studying the effect of an hydrophilic substrate 
since the surface of Vycor pores is well characterized^. 
In particular experimental studies with quasi-elastic neu- 
tron scattering and neutron resonance spin-echo^ on 
confined water indicated a slowing down of the dynamics 
with respect to the bulk and a study focused on the low 
hydration regime has evidenced upon supercooling the 
existence of a low frequency scattering excess typical of 
strong glass former—. 

In computer simulations of water confined in a pore 
of Vycor glass we found that strong layering effects are 
present where the molecules close to Vycor glass show 
very slow dynamics even at ambient temperatureiSiii*^^. 
In a more focused study on the low hydration regimei^ii^ 
we recently performed also a preliminary study of the 
residence time (RT) of the water molecules. We found 
that the RT is strongly dependent on the distance from 
the substrate and its distribution shows an anomalous 
non-brownian behaviour when the contribution of the 
molecules close to the substrate alone is considered. 

The aim of this paper is to show the connection be- 
tween the anomalous behaviour of the RT and the long 
time limit of the mean square displacement (MSD) in 
the low hydration regime. The paper is structured as fol- 



lows: in the next section we briefiy describe the system, 
give details of the computer simulation and discuss some 
structural properties of confined water— useful for the 
characterization of the RT behaviour. In the third sec- 
tion by considering the residence time distribution and 
the molecular diffusion we show the onset of an anoma- 
lous behaviour connected to the presence of the solid dis- 
ordered substrate. The last section is devoted to conclu- 
sions. 



II. STRUCTURAL PROPERTIES OF 
CONFINED WATER 

The molecular dynamics (MD) calculations have been 
performed in a cell of silica glass previously obtained by 
the usual quenching procedureii. 

Inside the glass cubic cell of 71 A a cylindrical cavity of 
40 A diameter has been carved. The surface of the cavity 
has been treated to reproduce the average properties of 
the surface of the pores of Vycor. Along this line hydro- 
gen ions (acidic hydrogens) have been added to oxygens 
non saturated by silicons in order to mimic the procedure 
followed by experimentalists in the preparation of the Vy- 
cor sample before hydration. The surface of silica glass 
can be considered as a prototype model for a disordered 
hydrophilic substrate. Being primarily interested in the 
dynamics of water the substrate is kept rigid. The water 
inserted in the cavity is simulated by using the SPC/E 
potential, where each molecule is represented by three 
charged sites. These sites interact also with the silicon 
and oxygen atoms of the substrate by means of an hy- 
drophilic potential described in previous workiii^. The 
molecular dynamics (MD) is performed at different hy- 
drations by varying the number of water molecules con- 
tained in the pore. For each hydration the system is 
equilibrated at different temperatures. The quantities of 
interest presented in the following are averaged over runs 
that extend up to 1.2 — 1.3 ns. We note that we have 
extended the simulation length at all temperatures with 
respect to ref— to improve statistics especially on the 
computation of the RT. 

The number of water molecules considered in this work 
are Nw = 500, Nw = 1000, Nw = 1500, which corre- 
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FIG. 1: For iV^ = 1500 water molecules at T = 300 K along 
the cylindrical radius R are shown: at the bottom density pro- 
files along the cylindrical radius R a.t T = 300 K normalized 
to the density of bulk water at ambient conditions; in the 
middle number of hydrogen bonds (HB) per molecule, the 
bold line represents the water- water HB, the dotted line the 
water- Vycor HB; at the top residence time along the cylindri- 
cal radius R. 
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FIG. 2; For Nw = 1000 water molecules at T = 300 K along 
the cylindrical radius R are shown: at the bottom density pro- 
files along the cylindrical radius _R at T = 300 K normalized 
to the density of bulk water at ambient conditions; in the 
middle number of hydrogen bonds (HB) per molecule, the 
bold line represents the water- water HB, the dotted line the 
water- Vycor HB; at the top residence time along the cylindri- 
cal radius R. 
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FIG. 3: For Nw = 500 water molecules at T = 300 K along 
the cylindrical radius R are shown: at the bottom density pro- 
files along the cylindrical radius i? at T = 300 K normalized 
to the density of bulk water at ambient conditions; in the 
middle number of hydrogen bonds (HB) per molecule, the 
bold line represents the water- water HB, the dotted line the 
water- Vycor HB; at the top residence time along the cylindri- 
cal radius R. 



spond to hydration levels of the pore of 19%, 38% and 
56% respectively, since the density corresponding to the 
full hydration in the experiments p — 0.878 g/cm^ is ob- 
tained in our geometry for Nw = 2600. 

The effects of the hydrophilic interaction of the sub- 
strate on the water molecules are shown in the bottom 
panels of Fig. ^ [3 and O for the different hydrations. 
The radial density profiles normalized to the density of 
bulk water at ambient conditions show the formation 
of two layers of molecules close to the substrate. The 
positions of the peaks of the double layer structure do 
not change with hydrations, the first layer is located at 
around R = 17.5 A and the second layer is at i? w 15.5 A 
with a minimum in the density at i? « 16 A. The heights 
of the peaks increase with increasing hydrations and for 
N]Y = 1500 the normalized density profile for confined 
water in the layers reaches values higher than one at am- 
bient temperature. 

In the middle panels of Fig. ^ [3 and |21 it is shown 
that the intermolecular hydrogen bond (HB) profiles in- 
crease and reach a maximum value at the position of 
the minimum of density {R « 16 A) where they start 
to go down in correspondence with the increase of the 
HB of water molecules with the atoms of the Vycor sur- 
face. The layering effect shown in the bottom panels is 
due to the formation of the Vycor- water HB. We found 
that the temperature has little effects on the density and 
HB profiles at all the hydration levels and for this reason 
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we show here only the results corresponding to ambient 
temperature. 

III. RESIDENCE TIME AND ANOMALOUS 
DIFFUSION 

At the top of Fig.lIllllandOlthe residence times (RT) of 
the water molecules at ambient temperature are reported 
along the pore radiusi^. The large oscillations of the RT 
appear modulated by the structure of the density profiles, 
reported in the bottom panels of the same figures. Apart 
for the molecules attached to the surface water resides 
for the longer time inside the shells, where the density 
profile reaches the highest values. The minima of the RT 
are located close to the minima of the density profiles. 

In Fig ^ are reported the calculations of the mean 
square displacement (MSD) at Nw = 1000 for decreasing 
temperatures. It is clear that after the ballistic regime at 
short time, at around 0.1 ps there is the onset of a cage 
effect characterized by the presence of a plateau which 
increases by lowering the temperature. The plateau is 
determined by the transient caging of the nearest neigh- 
bours. At longer times the MSD does not appear to 
reach the usual Brownian diffusion since the behaviour 
is sub-linear. Analogous results are found for the other 
hydration levels investigated. At this point further analy- 
sis is needed in order to clarify whether this subdiffusive 
behaviour is just a transient leading to normal brown- 
ian diffusion for longer times, although unreachable with 
normal computers, or it can be framed in the context of 
anomalous subdiffusive phenomena. 

Anomalous diffusion is generally speaking defined 
trough the time dependence of the MSD, the second mo- 
ment of the spatial coordinates of the diffusive particles, 
which generally has a long time dependence of the form 

<r2(i)>=ar (1) 

Anomalous diffusion corresponds to a ^ 1. In particular 
a > 1 is termed superdiffusion and a < 1 subdiffusion. 
From a theoretical point of view the origin of anoma- 
lous diffusion can be traced back to the analytic form 
of the distribution of the waiting times il^{t). Under the 
assumption that the amplitude of the random jumps is 
constant and finite anomalous diffusion has been shown 
to be generated by an inverse power law distribution for 
large timesiSiii 

iP{t) = At-f (2) 

In the case of ordinary brownian motion the long time 
limit of the distribution would decay with an exponential 
law. 

In our case the sublinear diffusion observed can be con- 
nected to the processes which take place close to the sub- 
strate and to the interaction of the water molecules with 
the disordered surfacei^. In this respect since oscillations 
of the RT appear so closely connected to the double layer 
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FIG. 4: MSD for Nw = 1000 at decreasing temperature: 
T = 300, 280, 260, 240, 220 K from above. 

structure it is of interest to look at the residence time dis- 
tribution of the water molecules close to the substrate. 

In Fig. O and El we report the residence time distribu- 
tions (RTD) at the highest T = 300 K and the lowest 
T — 240 K temperatures investigated for N\y = 1000 and 
Nw = 1500 . The RTD, calculated for the molecules in 
the layer 14 < i? < 20 A, shows the predicted power law 
behaviour of Eq.|21 while for the rest of the molecules we 
get an exponential decay, as shown in the inset of the fig- 
ures. The power law behaviour related to the temporal 
disorder of the distribution of the residence times of the 
molecules has been observed in computer simulation and 
experiments on water at contact with proteins^. They 
are specifically related to the interaction of the solvent 
with the protein sites. 

The power law decay of the RTD of the molecules in the 
6 A layer from the surface is determined by values of the 
exponents which are similar to the ones obtained for the 
RTD of water in few A shells close to protein hydration 
sites^. In particular for the case Nw = 1000 we have 
^ = 1.54±0.05 at T = 300 K and = 1.50±0.05 at T = 
240 K, while for Nw = 1500 the fits yield ^ = 1.50±0.05 
at T = 300 K and = 1.52 ±0.05 at T = 240 K. We note 
that the present result for /i at Nw — 1500 and ambient 
temperature is slightly different from the preliminary one 
reported in ref^, where the statistics was poorer. 

The sublinear behaviour of the MSD is connected to 
the power law decay of the RTD by the asymptotic tem- 
poral dependence 

< r2(t) >oct^-i (3) 

From the fit of the long time behaviour of the MSD 
for Nw = 1000 and Nw = 1500, reported in Fig. □ and 
Fig. |S1 we get the values of the exponent fi which are 
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FIG. 5: Log-log plot of the residence time distribution (RTD) 
of the water molecules in the layer 14 < 7? < 20 A for Nw = 
1000 at temperatures T — 300 K (open triangles) and T = 
240 K (open squares). The fits (bold lines) are done with a 
power law At'". A = 0.059, fi = 1.54 ± 0.05 at T = 300 K. 
A = 0.021, n = 1.50 ± 0.05 at T = 240 K. In the inset are 
shown in a linear-log scale the RTD for the molecules in the 
layer < J? < 14 A. 
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FIG. 7: MSD of water molecules in the layer 14 < i? < 20 A 
for Nw = 1000 at temperatures T = 300 K and T = 240 K 
from above. The long dashed lines are the fits to a sublinear 
behaviour < >(x t" with a = 0.46 ± 0.05 at T = 300 K 
and a = 0.48 ± 0.05 at T = 240 K. In the inset are reported 
the functions <r'^ > jt. 
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FIG. 6: Log-log plot of the residence time distribution (RTD) 
of the water molecules in the layer 14 < J? < 20 A for Nw = 
1500 at temperatures T — 300 K (open triangles) and T = 
240 K (open square). The fits (bold lines) are done with a 
power law At-t". A = 0.059, fi = 1.50 ± 0.05 at T = 300 K. 
A = 0.026, II = 1.52 ± 0.05 at T = 240 K. In the inset are 
shown in a linear-log scale the RTD for the molecules in the 
layer < 7? < 14 A. 



consistent with the exponents obtained from the power 
law behaviour of the RTD seen in Fig.|31and Fig. El 

In the lower hydration case Nw = 500, Fig. El the 
RTD decays for long time with an exponent similar to 
the previous cases, fi — 1.45 ± 0.05 and ^ = 1.55 ± 0.05 
for T = 300 K and T = 240 K respectively, but the MSD 
show, as seen in the inset of Fig. El a long time behaviour 
not in agreement with the one predicted by Eq. El 

The asymptotic behaviour of the MSD shows a further 
slowing down with respect to the higher hydrations. This 
behaviour is likely connected to the fact that the water 
molecules are arranged in clusters close to the substrate. 



IV. CONCLUSIONS 

The dynamical properties of confined water are ex- 
pected to be changed by the interaction with the sub- 
strate. We performed computer simulation of water 
molecules confined in a silica pore in low hydration 
regimes, where the larger amount of water resides in the 
shells closer to the hydrophilic surface. We found drastic 
changes of the diffusion of the water molecules. From a 
layer analysis for the investigated hydrations we find that 
the diffusive regime at long time of the molecules close 
to the substrate is characterized by a sublinear trend. 

Anomalous diffusion phenomena are related to a tem- 
poral disorder typical of particles which diffuse close to 
and interact with a disordered surface. Different inter- 
action processes between the water molecules and the 
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sites of the substrate modulate the residence time of the 
molecules^. A dispersive transport regime related to tem- 
poral disorder shows up in the power law decay of the 
residence time distribution with an exponent which de- 
termines also the long time tail of the mean square dis- 
placement. In our system the exponent of the long time 
behaviour of the mean square displacement is related to 
the long time decay of the residence time distribution of 
the molecules in the same layer for the cases Nw = 1500 
and = 1000 as theoretically predicted. For the lowest 
hydration case {Nw = 500) the mobility of the molecules 
is more strongly modulated by the substrate with respect 
to the higher hydrations. The formation of clusters of 
molecules close to the solid surface does not appear to 
modify the long time decay of the RTD but it induces a 
further slowing down of the dynamics with a violation of 
the expected behaviour of the MSD. 



FIG. 8: MSD of water molecules in the layer 14 < i? < 20 A 
for Nw = 1500 at temperatures T = 300 K and T = 240 K 
from above. The long dashed lines are the fits to a sublinear 
behaviour < >oc t" with a = 0.45 ± 0.05 at T = 300 K 
and a = 0.48 ± 0.05 at T = 240 K. In the inset are reported 
the functions < > /t. 
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FIG. 9: Log-log plot of the residence time distribution (RTD) 
of the water molecules in the layer 14 < 7? < 20 A for Nw = 
500 at temperatures T — 300 K (open squares) and T — 240 K 
(open triangles). The fit (bold lines) is done with a power law 
At-'', with A = 0.048, fj. = 1.45 ± 0.05 for T = 300 K and 
A = 0.041, ^ = 1.55±0.05 for T = 240 K. In the inset MSD of 
water molecules in the layer 14 < ii < 20 A for A^vi^ = 500 at 
temperatures T = 300 K and T = 240 K. The bold hue is the 
fit to a sublinear behaviour < >(x t"' with a — 0.34 ± 0.05 
for both T = 300 K and T = 240 K. 
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